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The Kondo effect, one of the oldest correlation phenomena known in con- 
densed matter physics [Ij, has regained attention due to scanning tunnel- 
ing spectroscopy (STS) experiments performed on single magnetic impurities 
[21 |3]. Despite the sub-nanometer resolution capability of local probe tech- 
niques one of the fundamental aspects of Kondo physics, its spatial extension, 
is still subject to discussion. Up to now all STS studies on single adsorbed 
atoms have shown that observable Kondo features rapidly vanish with increas- 
ing distance from the impurity |4|j [H |6l [71 El |9] . Here we report on a hitherto 
unobserved long range Kondo signature for single magnetic atoms of Fe and 
Co buried under a Cu(lOO) surface. We present a theoretical interpretation 
of the measured signatures using a combined approach of band structure and 
many-body numerical renormalization group (NRG) calculations. These are 
in excellent agreement with the rich spatially and spectroscopically resolved 
experimental data. 
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The interaction of a single magnetic impurity with the surrounding electron gas of a 
non-magnetic metal leads to fascinating phenomena in the low temperature limit, which 
are summarized by the term Kondo effect [1]. Such an impurity has a localized spin 
moment that interacts with the electrons of the conduction band. If the system is cooled 
below a characteristic temperature, the Kondo temperature T^, a correlated electronic 
state develops and the impurity spin is screened. The most prominent fingerprint of this 
many body singlet state is a narrow resonance at the Fermi energy 6:i? in the single particle 
spectrum of the impurity, called Kondo or Abrikosov-Suhl resonance. The existence of 
this Kondo resonance has been experimentally confirmed for dense systems with high 
resolution photoemission electron spectroscopy and inverse photoemission [TOl [H] . Due 
to their limited spatial resolution these measurements always probe a very large ensemble 
of magnetic atoms. With its capability to study local electronic properties with high 
spatial and energetic resolution. Scanning Tunneling Spectroscopy (STS) has paved the 
way to access individual impurities [21 [3] . 

A theoretical prediction for the local density of states (LDOS) - the key quantity 
measured in STS experiments - was first provided by Ujsaghy et al [12j. According to 
their calculations the Kondo resonance induces strong spectroscopic signatures at the 
Fermi energy whose line shape is oscillatory with distance to the impurity. Since the first 
STS studies in 1998 0|3] a lot of experiments on magnetic atoms and molecules on metal 
surfaces were performed, all revealing Kondo fingerprints [51 |6l [71 El |9] . However, it is 
worth noting that all previous STS experiments on isolated ad atoms have reported that 
the Kondo signature rapidly vanishes within a few angstrom and no variation of the line 
shape occurs (for a review on ad atom Kondo systems see [I3]). 

In the present work we follow a novel route and investigate single isolated magnetic 
impurities buried below the surface with a low temperature STM operating at 6K. It has 
been recently shown that the anisotropy of the copper Fermi surface leads to a strong 
directional propagation of quasi particles, which is called electron focusing [14j. This 
eflFect gives access to individual bulk impurities in a metal that were previously assumed 
to be "invisible" due to charge screening. Following these observations, dilute magnetic 
alloys were prepared on a clean Cu(lOO) single crystal by adding a small amount (0.02%) 
of either Co or Fe to the topmost monolayers. We have chosen Co and Fe because of their 
different Kondo temperature. This allows to test the universal character of the Kondo 
effect. 
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Figure 1: Variation of the LDOS with the lateral distance of the tip. a) STM topography 
(1.8nm X 1.8nm, 2nA) of a 4th layer Fe impurity for different bias voltages V. For 

V < the LDOS directly above the impurity is reduced. By increasing V the central ring 
contracts to a plateau like maximum and the size of the pattern decreases, b) Section 
of AdI/dV = dl/dV — dIo/dV along the [010] direction. The tip height was adjusted to 
give a dIo/dV of 13.3pA/mV for the free surface. An energy dependent phase shift of the 
interference pattern can be observed. The vertical black line would resemble the phase 
front for a energy independent scattering phase and indicates that the overall phase shift 
caused by the resonance is less than tt. Single spectra (purple curve) for four different 
lateral distances (marked in b by black arrows) are depicted in c). To illustrate the relative 
amplitude of the Kondo signal, the data is normalized to the differential-conductance of 
the free surface (blue scale on the right side). Black curve: Calculated spectra obtained 
by fitting a phenomenological expression of the Kondo resonance to the dl/dV signal (see 
supplementary information) . 

As a first striking example how the Kondo effect influences the energy dependent 
scattering behavior on the mV scale, figure la shows STM topographies of a fourth layer Fe 
impurity for different bias voltages near zero mV. The local minimum of the LDOS present 
in the center of interference pattern ioi V < develops into a plateau like maximum for 

V > 0. In figure lb) the differential conductance as a function of bias voltage and one 
spatial coordinate y across the impurity pattern is shown. The crossover observed in 
the topographies occurs very close to zero bias. In figure Ic) four spectra for different 
positions A-D are shown, illustrating very clearly that a single Kondo atom buried below 
the surface of copper induces long range spectral signatures that depend on the distance 
to the impurity. A second possibility to investigate the Kondo effect versus distance is 
to look at impurities situated in different depth d below the surface. Figure 2 shows as 
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an example STS-data of Co atoms, comparing single spectra (purple curves) measured 
directly above the impurities = 0). The lateral variations of dl/dV are again depicted 
as sections (upper part of figure 2). All reveal a constriction of the pattern for positive 
bias voltages. The comparison between Co and Fe data demonstrates that both impurity 
species show similar behavior on completely different energy scales (for instance compare 
the 4th layer Fe in figure lb with a 4th layer Co in figure 2). 




Vbltage[mV] Vbltage[mV] U)ltage[mV] U)ltage[my] Vbltage[mV| 



Figure 2: STS-data of subsurface Co impurities in 3 to 7 layers below the Cu(lOO) sur- 
face. Spectrum sections AdI/dV = dl/dV — dIo/dV along the [010] direction are are 
shown in the upper part. A strong change in the scattering behavior can be seen for 
all impurity depths. The black bar corresponds to a length of 0.5nm for every section. 
Single AdI/dV spectra (purple curve) measured with a tip position direct above the im- 
purity {d = 3..7ML) are depicted in the lower part. The tip height was adjusted to give 
a differential-conductance for the free surface dIo/dV of 3.2pA/mV (3ML, 4ML, 6ML), 
5.3pA/mV (5ML) and 6.4pA/mV (7ML). To illustrate the relative amplitude of the Kondo 
signal, the data is normalized to the differential-conductance of the free surface (blue scale 
on the right side). Black curve: Calculated spectra obtained by fitting a phenomenological 
expression of the Kondo resonance to the dl/dV signal (see supplementary information). 

For a quantitative analysis, we have to use advanced tools of quantum many-particle 
theory as the Kondo problem is a genuine many-body effect. The well-known universal 
behavior of the Kondo effect and its associated low-energy fingerprints allows to apply the 
Single Impurity Anderson Model. Within this model the localized orbital of the impurity 
is described by (i) a single level that couples to non-interacting conduction band electrons 
and (ii) a Coulomb interaction U between electrons on the site. 

The effect of the Kondo resonance in the spectral function can be understood in 
analogy to other fields of physics, e.g. scattering of electrons at a potential well: the 
resonance causes an enhanced scattering amplitude and an energy dependent phase shift 
of (quasi-) particles near the resonance energy. 
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In many-body formalism the measured LDOS is the imaginary part of the single 
electron Green's function G(x, x, s). To calculate this quantity we use the Dyson equation, 
which connects the Green's function and hence the LDOS of the perturbed system to the 
Green's function G^(x, x', 6:) of the unperturbed conduction band electrons via a T-Matrix 
T (ximp, s). Approximating the impurity to be a point scatterer at position Xi^p the LDOS 
change is given by 

ALDOS(x,e) = --Im [^^(x, Xi^p, e)r (xi^p, e) ^^(xi^p, x, e)] (1) 

TT 

Calculating the band structure of copper using a linear combination of atomic orbitals 
(LCAO) approach (figure 3a) and treating the surface as a potential step by fixing the con- 
tinuity conditions for the wave functions we obtain the free Green's functions G^(x, x^, s) 
for copper (figure 3b and supplementary information). For the evaluation of the T-Matrix 
we apply the numerical renormalization group (NRG) [15j. At this point the imaginary 
part of the unperturbed Green's functions G°(ximp, Xi^p, 6^) at the impurity position is 
needed to describe the states of the conduction band that couple to the impurity. Note 
that the presence of the surface leads to oscillations in this quantity (see figure 3c) , which 
can be understood as a standing wave pattern of electrons being refiected at the surface. 
The NRG calculations result in a Kondo resonance, in the T-Matrix near Sjp (figure 3d) 
with a resonance height and width showing similar oscillatory behavior as seen in the 
unperturbed Green's functions G°(x, x',6:). 

With the T-Matrix obtained by the above procedure we calculate the LDOS change 
using equation (1). Values for energy and hybridization of the localized orbital were taken 
from ab-initio results [16j. The only free parameter left, the Coulomb interaction [/, was 
adjusted by comparison to the experiment. The calculated LDOS sections (figure 3e) are 
in excellent agreement with the measured data (compare figure 2). The experimentally 
observed periodicity (compare a 3rd and 7th layer impurity) is also found in the simulation. 

The rich spatial and spectroscopic information of the measured interference patterns 
allow further investigation of the Kondo resonance. The resonance width is proportional 
to the Kondo temperature (see supplementary information). We fit the experimental data 
to a phenomenological form found by Frota (see discussion further down) for all lateral 
positions and a constant impurity depth. Averaging over all different depths we get a 
Kondo temperature Tk for Fe impurities of 32(16)K and for Co of 655(155)K. Results 
from macroscopic bulk measurements ^17j give similar values. 
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Figure 3: Road map of the theoretical modeh a) band structure of copper calculated 
using LCAO and b) the related free propagator, c) Unperturbed LDOS as function of 
depth including effects of the surface, d) The Kondo resonance for impurities located 
at different positions, e) Calculated ALDOS section along the [010] direction for Co 
impurities situated in 3 - 7 monolayers below the surface. Single spectra taken directly 
above the impurity position are shown below. 



Investigating the Kondo temperature as a function of the lateral position for a constant 
impurity depth we found an unexpected variation of the resonance width. As an example 
the 4th layer Fe impurity (figure Ic) shows a variation from Tk (A) = 29(9)K, Tk (B) = 
35(24)K, Tk (C) = 42(18)K and Tk (D) = 44(26)K from position A to D. This behavior 
is not included in our theory so far. One possible explanation employs the actual orbital 
structure of the d-level in Fe or Co, which leads to a dependence of the T-Matrix on the 
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propagation direction. 

Going further, we investigate not only the long range signature of the scattering am- 
plitude but also its phase. During our analysis it turned out that the widely used (com- 
plex) Lorentzian approximation of the Kondo resonance does not give the best descrip- 
tion to both our experimental and theoretical data (see supplementary information). A 
Lorentzian for the T-Matrix results in the often used Fano line shapes [18]. Better fits 
were obtained by a phenomenological form found by Frota [191 EQ] • The scattering ampli- 
tude of the Kondo resonance decays much weaker with energy than a Lorentzian, which 
can be viewed as a minor correction. More significant, the shift of the scattering phase 
due to the Kondo resonance is strongly overestimated as a Lorentzian always results in 
a phase shift of tt. In the range of energies considered here, our NRG calculations give a 
phase shift of only 7r/2 which is also descried by the phenomenological form of Frota. As 
one can see in figure lb the phase shift measured in experiment is also smaller than tt. 

The mapping of the scattering phase as a function of energy opens a new way to 
discriminate between single particle resonances and Kondo scattering. Our approach 
allows further investigations of the properties of magnetic impurities, like the splitting 
of the resonance by an applied magnetic field. This was recently shown on adsorbed 
atoms [21] and should also effect the scattering behavior of subsurface impurities. With 
available magnetic fields manganese may be a good candidate for such an experiment. 
Since we observe a long range Kondo signature this opens a new way to examine the 
interaction between two or more Kondo atoms with each other, or the effect of an interface 
in real space. Finally this observation may give an access to one of the most controversial 
discussed terms in Kondo physics: the meaning and size of the "Kondo cloud". With the 
long range LDOS signatures of buried impurities charge density oscillations are accessible 
and the term " Kondo cloud" may be defined in a way which is observable and consistent 
with experiments. 
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1 Supplementary Information: Long range Kondo 
signature of a single magnetic impurity 

1.1 Sample preparation 

The Cu(lOO) single crystal substrate is cleaned by several cycles of argon bombardment 
and heating. On this substrate thin films were prepared by simultaneous deposition of 
copper and the impurity material from two electron beam evaporators, leading to copper 
alloys with a small amount (0.02%) of iron or cobalt embedded. The measurements were 
performed with a home built low temperature STM operating in ultra high vacuum with 
a base pressure better than 5 • 10~^^ mbar. We use electrochemically etched tungsten tips, 
prepared by annealing and argon ion bombardment. The performance of the tip is tuned 
by controlled voltage pulses and smooth sample tip contacts. All standing wave patterns 
are ordered to the lateral size. The depth of the impurity is a monotone increasing function 
of the size of the standing wave pattern [H]. resolved topographies, which exists for most 
of the data presented in this work (for example see figure SI). Since the Cu(lOO) has a 
fee crystal structure and Fe and Co atoms are situated on substitutional lattice sites, the 
impurity contrast has to follow a certain ordering with respect to the host lattice. The 
standing wave pattern of an odd layer impurity is centered on an corrugation maximum 
of the surface. If the impurity is positioned in an even layer the center is located between 
the corrugation maxima. 

1.2 STS data acquisition 

To investigate the Kondo eflFect and its infiuence on the LDOS (?/, s) we use STS spectro- 
scopic data acquired by recording an I(U) curve at every scanning point with interrupted 
feedback loop. Further data processing included averaging and numerical diflFerentiation. 
This provides a complete map of the differential conductance dl/dV as a function of lat- 
eral tip position and bias voltage on the Cu(lOO) surface. Since the tip is stabilized in 
constant current mode, the STS data is affected by different tip-sample distances. This 
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artifact can be successfully removed by normalizing the STS data to a constant sample 
tip distance [22j. To extract the LDOS change we subtract the differential conductance 
of the free surface dl/dV. This quantity is proportional to the change in LDOS j23] at 
the lateral tip position. In order to perform a quantitative analysis of the experimental 
data, drift effects which occur during acquiring the spectroscopic data, were removed by 
comparing the simultaneously measured topography with a calibrated one. 

1.3 Equation of motion technique for the SI AM 

For the simulation of the studied system we have chosen the single impurity Anderson 
model (SIAM) [24]. Within the SIAM the impurity is modeled by an effective single level 
with energy Cd- For double occupancy the electrons are exposed to the repulsive coulomb 
interaction U. The hybridization of the impurity to the copper crystal is given via the 
hybridization parameter Vk that describes the probability of the hopping of an electron 
from a conduction electron state with momentum k to the impurity level e^. Altogether 
the Hamiltonian for the SIAM in momentum space reads 

^^SIAM = V //^^^ + //'^P 
ka 
ka 
a 

where and c^^ denote the creation and annihilation operators for an electron with 
momentum k and spin a = {t, i} ( d denotes the d(f)-state of the impurity), n^^ = c\^c^^ 
the number operator and e{k) the dispersion relation of the free system. Setting up the 
equations of motion, the Green's function Gkak'a'{^) can be written via the T-Matrix [Ij 
(G^^;./^/ are diagonal in k and a) 

G,k'{e) = GUe) + E GUe)TUe)GlAe) 

Im 

with the free Green's function G^^j^,^,{e) = -^^^^^f^- For the SIAM the T-matrix is given 
by: 

TkA^) = V:Gdd{e)Vk' (2) 
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Transforming this equation to real space yields: 



G(x, x', e 



:) = G" 



(x, x', e 



) + 



// 



(x,z,e 



)T(z,z',e 



:)G°(z',x',e 



) 



(3) 



with 



(4) 



In the foUowing we focus on a true point-scatterer at position x^. 



leading to 



G(x, x', e) = G^(x, x', e) + G^(x, x^^^, e)V'^Gdd{^^ ^imp)G^{^imp, x', e) (5) 



In the last equation we have inserted ^imp into the arguments of Gdd{^^ ^imp) to underline 
that the system is not truly decoupled as Odd still depends on the free Green's function 
G^i'^imp^ ^imp) ^) ^t the impurity position. The free Green's functions are calculated via a 
Linear Combination of Atomic Orbitals (LCAO) method including the surface and Gdd{^) 
is calculated by applying the numerical renormalization group (NRG). 

1.4 Evaluating the free Green's function via LCAO 

The Propagator G^{x^x'^e) of the unperturbed conduction band electrons was obtained 
using spectral representation from the band structure E\^ i, of copper and the wave func- 



This quantity describes the propagation of a single electron with energy e from a point 
source at x' to other positions x in the system. The band structure E\^ i, of copper was 
calculated by an LCAO approach with parameters taken from [25j. The wave functions 
^k,z/ were derived from plane waves while reflections of the electrons by the surface and 
a realistic decay into the vacuum were taken into account. This decay depends on the 
energy and the wave vector component parallel to the surface: 



Here $ = A.GbeV is the work function and the energy E is measured with respect to the 
Fermi-Level. The surface is orientated perpendicular to the z-axis and deflnes the border 



tions ^^kz/^ 






(7) 
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between the two adjacent subspaces crystal {z <) and vacuum (z > 0). As the fee band 
structure is invariant under reflections at the (OOl)-plane, the wave functions are: 

^^k,.(x||, z) = (1 + M'y'^' e^*^""" [e''-' + rk,.e-^^^^] z<0 (8) 

^k,.(x||, = (1 + |rk,,|')"'^' 6^*^11^11 (1 + rk,.) e-^^-' z>0 (9) 

with the reflection coeflicient rk,zy = {ik^ + i^i^^u) / {ik^ — /^k,?/) providing continuity of the 
wave-functions and their spatial derivative at the surface. The LDOS of the unperturbed 
system in both subspaces is then: 



go{z,s) = y^ / d^k5{s- E^^^)\^^^^{^\\,z)f 

„ Jl.BZ 



(10) 



The STM-Tip always probes the LDOS in the vacuum and the impurity is always lo- 
cated within the crystal. Therefore only the Green's function connecting a position 
= (0,0, —rf) inside the crystal with a position x = in the vacuum is of in- 

terest. Here d is the depths of the impurity below the surface and h is the the tip to 
sample distance. This particular Green's function is given by: 

/ d'k-^^^^^^^ (11) 

, Jl.BZ (l + K.|')(£-^k,. + ^0+) ^ ' 

The Propagator for the opposite direction is identical due to time reversal symmetry. 
Hence the vacuum LDOS at a distance h above the surface which is modifled by an 
impurity in depth d below the surface is given by: 

p^(x||,/i,s) = p^ih,e) - -Im [G'{^\\,h,d,e)'T{e)] (12) 

The LDOS 6:) of the unperturbed crystal enters the NRG calculations in order to 



calculate the T-Matrix for the Dyson equation (12). In the main article (see flgure 3b) the 
vacuum LDOS p^(x||, /i, s) of the perturbed system is shown for h = 7A and d = 3...7ML. 

1.5 Evaluating the Green's function of the impurity with the 
NRG 

For the evaluation of the Green's function of the impurity Gddi^i^imp) we apply the 
Numerical Renormalization Group (NRG). The SIAM can be interpreted as an impurity 
coupled to a bath of non-interacting states. The coupling is sufliciently described by the 
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hybridization function A(6), that is proportional to the density of conduction states that 
hybridize with the impurity state. For a point-hke hybridization in real space this is given 
by the LDOS at the impurity position x^^^ and the hopping probability V. Therefore the 
hybridization function is given by: 

A(e) = 7rV^LB0S{^,mp) = -V^3mG'{^,mp. ^,mp. e) (13) 

The free Green's function at the impurity position is given by the LCAO calculations. 
Within the NRG the hybridization function is logarithmically discretized and mapped on 
a semi-infinite chain with the impurity at one end. The logarithmic discretization ensures 
that the hopping amplitudes of the chain fall off exponentially. This allows an iterative 
diagonalization of the chain, whereby at every step one site is added to the chain. As 
the increasing Hilbert space of this procedure does not allow a complete diagonalization 
one has to introduce a truncation scheme in the spirit of the renormalization group. In 
the NRG only the A^^ lowest lying eigenstates are retained and used to build up the 
next Hamiltonian, thus keeping the size of the Hilbert space constant. The result of the 
iterative diagonalization scheme are the many-particle energies En{r) with r = 1, A^. 
For the evaluation of the Green's function we use the "self-energy trick" that maintains 
the Friedel sum rule for the peak height of the Kondo resonance (for a review see 
The SIAM has three free parameters U and V. From [12l [16] we take Cd = —O.SeV 
and A(ei?) = 0.2eV^ so that hybridization is extracted accordingly to he V = 0.447. The 
interaction strength U is adjusted to the experimentally observation. The best accordance 
for cobalt is given by a value oi U = 1.2eV and for iron for U = 2AeV. On should point 
out that for cobalt we are not in the true Kondo regime that is riven . > 3. The 
resulting spectral functions for the impurity show the Kondo resonance close the Fermi 
energy. The width and height of the resonance show similar oscillatory behavior as the 
hybridization function. The oscillations in the peak height are in accordance with the 
Friedel sum rule Gdd{^F,^imp) = ^A(e/,x,^p) • 

1.6 Extracting the Kondo temperature 

In this work we extract the Kondo temperature from the half-width at half maximum 
(HWHM) of the Kondo resonance. We use Wilson's 1975 definition [26j of the Kondo 
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temperature T^w- 

Tkw = D{Jp)exp 



In \2Jp\ - 1.5824|2 Jp| + 0( JV) 



(14) 



2Jp 2 

^ DV2\Jp\exp{-l/\2Jp\) (15) 

Extracting the Kondo Temperatures from the Kondo resonance is done via the HWHM 
In the crudest approximation the Kondo resonance is modeled via an Lorentzian 
Pl{^) oc ff^- Actually this is only a valid approximation in the asymptotic e ^ 
region [271 [28j . The HWHM is then directly given via the = F^^. A better line shape 
gives a phenomenological form found by Frota et al.[20], [I9j. A comparison between the 
phenomenological function found by Frota, the conventionally used Lorentzian and the 
resonance calculated with the NRG are shown in figure [4} One may clearly see the better 
modeling of the resonance by the Frota form. 



PF{e) = Im 



1 iF, 



i V ^Fi^ + £ 



(16) 



The HWHM of this form is given by A^? = 2.542 • F^?. Therefore setting this form on the 
same value at the HWHM as the Kondo resonance gives = 0.393 • A^. Following [27 \ 
HWHM is proportional to the Kondo resonance. 

^z^tko ^ 3.7 'Tkw (17) 

Therefore the extracted parameter = 0.393 -3. 7 kBTxw = 1.455/c^Txm/ is taken for the 
Frota form. Using this proportionality constant, one has to remind that it is only valid in 
the wide band limes for a constant band. For the regarded system this is not the case as 
the presence of the surface gives small deviations. Therefore the published values for the 
Kondo Temperatures for the mean values of the resonance widths are to be interpreted 
with caution. 



1.7 Phenomenological expression for the LDOS 

To extract the key quantities from the experiment and to focus on the universal features 
of the Kondo resonance in the low energy regime, we use the former mentioned phe- 
nomenological form for the T-Matrix. The free propagation of band electrons is modeled 
by an energy-independent phase factor G^{r) ^ 6^^9(0/2 ^j^]^ ^ ^ j-^^ _ Xjj^pl. The energy 
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Figure 4: Comparison of the Kondo resonance calculated with the phenomenological Frota 
form and Lorentzian. 

independence is justified as the energy is nearly constant around the Fermi energy (see 
figure 4c). The frequency of the oscillation in the free Green's function is given by the 
parameter (/>g(r) = 2/cf • r, which determines the line shape. Using these two approxima- 
tions we finally obtain the following fit function from equation ^ to describe the LDOS 
change due to an impurity in distance r. 



Here Sk is the position of the Kondo resonance and F^? is proportional to the resonance 
width (see former section). A linear voltage lope /56 + 7 is added to account for additional, 
approximately energy independent background scattering processes. 
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